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Abstract 

We evaluate the thermodynamic critical angular velocity Vlc{T) for creation of a 
vortex of lowest quantized angular momentum in a strictly two-dimensional Bose 
gas at temperature T, using a mean-field two-fluid model for the condensate and 
the thermal cloud. Our results show that (i) a Thomas-Fermi description of the 
condensate badly fails in predicting the particle density profiles and the energy of 
the vortex as functions of T; and (ii) an extrapolation of a simple Thomas-Fermi 
formula for r2c(0) is nevertheless approximately useful up to T ~ O.STc. 
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1 Introduction 



Quasi-two-dimensional (2D) Bose-Einstein-condensed gases (BEC) have been 
attracting increasing attention in recent years. Flatter and flatter condensates 
have been produced inside magneto-optical harmonic traps by squeezing the 
anisotropy parameter, measured as the ratio between the radial and axial 
trap frequencies (see Gorlitz et al. [1] and references therein). As the gas 
approaches the 2D hmit, its coUisional properties start to influence the boson- 
boson coupling parameter, which becomes dependent on the gas density [2,3]. 
Understanding the behaviour of vortices in this regime is important, since they 
reflect the superfluid nature of the condensate [4-6] and are expected to play 
a role in the transition from the superfluid to the normal state [7]. 

An experimental method for the creation of quantized vortices in a trapped 

BEC has made use of an "optical spoon" [8,9], whereby the condensate in an 
elongated trap is set into rotation by stirring with laser beams. This exper- 
iment on a conflned boson gas is conceptually the analogue of the rotating 
bucket experiment on bulk superfluid Helium, with the difference that the 
thermal cloud is inhomogeneously distributed and lies mostly outside the con- 
densate cloud. A quantized vortex flrst appears in the condensate at a critical 
angular frequency of stirring which corresponds to an instability of the vortex- 
free state. A lower bound for the critical frequency can be assessed from the 
energy of a vortex, defined as the difference in the internal energy of the gas 
with and without a vortex. Observation of a vortex in a trapped gas is difficult 
since the vortex core is small in comparison to the size of the boson cloud, 
but the size of the core increases during free expansion and indeed vortices 
were first observed experimentally by releasing the trap and allowing ballistic 
expansion of the cloud [8] . 

In this Letter we consider a 2D rotating condensate at finite temperature 
and evaluate the particle density profiles and the energy of a vortex within 
a strictly 2D model for the boson-boson coupling. This is appropriate to a 
situation in which the s-wave scattering length starts to exceed the vertical 
confinement length. Previous work has established that the dimensionality of 
the scattering collisions strongly affects the equilibrium density profiles [10,11] 
and the process of free expansion of a BEC containing a vortex [12]. A similar 
study of the density profiles of a rotating BEC in 3D geometry at finite tem- 
perature has been carried out by Mizushima et al. [13], who also determined 
the location of various dynamical instabilities within the Bogoliubov-Popov 
theory. 

We begin, therefore, by introducing our description of a 2D BEC containing 
a vortex at finite temperature T. This uses a mean- field two-fiuid model for 
the condensate and the thermal cloud. 
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2 The model 



The BEC is subject to an anisotropic harmonic confinement characterized 
by the radial trap frequency u!± and by the axial frequency Xu!± with A ^ 1. 
Motions along the z direction are suppressed and the condensate wave function 
is determined by a 2D equation of motion in the {x, y} plane. 

The order parameter for a 2D condensate accommodating a quantized vortex 

state of angular momentum Hk per particle is written as $(r) = i(j[r) exp^incj)), 
with the azimuthal angle. The wave function i/j{r) then obeys the nonlinear 
Schrodinger equation (NLSE) 



h!^ h^K^ 1 

"2^^'^' ^ 2^ ^ -m^ir^ + g^n^{r) + 2^inr(r) 



'^(r) = iJ>ip{r) (1) 



where fi is the chemical potential, m the atomic mass, nc(r) = |'?/'(r)p the con- 
densate density, and nT{r) the density distribution of the thermal cloud. The 
2D coupling parameters gj, with j — 2 for condensate-condensate repulsions 
and j — 1 for condensate-noncondensate repulsions, are given by 

\n\4hy{jmi^a^)\ ' ^ ^ 

where a is the 5- wave scattering length [2,3]. In Eq. (2) we have omitted a 
term due to thermal excitations [11], which is negligible in the temperature 
range of present interest (T < 0.5Tc, with T^. the critical temperature). 

In our model the atoms in the thermal cloud are not put directly into rotation, 

but feel the rotating condensate through the mean-field interactions. In the 
Hartree-Fock approximation [14] the thermal cloud is treated as an ideal gas 
subject to the effective potential 

^eff{r) = ^muilr^ + 2ginc{r) . (3) 
The density distribution of the thermal cloud is then given by 



TTX I 



(4) 



with 13 — l/{kBT). The momentum cut-off po in Eq. (4) is a simple expedient 

to eliminate discontinuities in the density profiles that can occur near the 
Thomas-Fermi radius of the condensate. We take po = 2y^m^7^ (see for 
instance [15]) which is equivalent to adding the term 2ginT to the effective 
potential in Eq. (3). 



3 



We solve self-consistently the coupled Eqs. (l)-(4) together with the condition 
that the areal integral of nc(r)+ny(r) is equal to the total number of parti- 
cles. The differential equation (1) is solved iteratively by discretization, using 
a two-step Crank-Nicholson scheme [16]. In Sec. 3 we shall also compare the 
results with those obtained in the Thomas-Fermi approximation by dropping 
the radial kinetic energy term in Eq. (1). 



2. 1 The energy of a vortex 



The critical angular velocity measured in experiments where vortex nucleation 
occurs from a dynamical instability [8] depends strongly on the shape of the 
perturbation. However, a lower bound for the angular velocity required to 
produce a single-vortex state can be estimated once the energies of the states 
with and without the vortex are known. Since the angular momentum per 
particle is hn, the critical (thermodynamic) angular velocity is given by 

This expression follows by equating the energy of the vortex state in the 
rotating frame, that is — ilcL^, to the energy i?K=o of the vortex-free state. 

In the noninteracting case at zero temperature, the energy difference per parti- 
cle is simply Hkuj^, so that Q^. is just the trap frequency in the {x, y} plane. For 
the interacting gas at zero temperature in the Thomas-Fermi approximation, 
Eq. (5) reduces to the expression [17,4] 

Here R = and ^ = Rfiu;±/2ijL are the Thomas- Fermi radius and 

the healing length respectively. 

In the general case of an interacting gas at finite temperature, we have to 
evaluate numerically the total energy as the sum of four terms [18], 

Ek = E\iiD,c + -E'trap + -E'int + -E'kin.T • (7) 

These terms are the kinetic energy of the condensate 

b.,.,c^/<<v^-m(-|1v^ + 1^1)^(0. (8) 

the energy of confinement 

^trap = ^mul J cPrr'^[nc{r) + ririr)] , (9) 
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the interaction energy 

^int = ^ / d^r[g2nl{r) + 4^in,(r)nr(r) + 2gin^^{r)] , (10) 
and the kinetic energy of the thermal cloud 

£^kin,T = Jd"^J^2^ {exp[/3(pV2m + ^//(r) - n)] - l}"' . (11) 

In Sec. 3 we compare the angular frequency obtained by the calculation of 
the total energy of the gas with and without a vortex with those obtained 
from a Thomas-Fermi calculation and from extrapolating Eq. (6) at finite 
temperature through the temperature dependence of the chemical potential. 



3 Results and discussion 

For a numerical illustration we have taken k — 1 and chosen values of the 
system parameters as appropriate to the gas of ^^Na atoms studied in the 
experiments of Gorlitz et al. [1], namely u± — 2x 188.4 Hz, a = 2.8 nm, and 
= 5 X lO''. Wc arc implicitly assuming, however, that the trap has been 
axially squeezed to reach the strictly 2D scattering regime. 
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Fig. 1. Density profiles n(r) for the condensate and the thermal cloud (in units of 

a^/N, with afio = ^Jhjmiij\)j versus radial distance r (in units of aho) at tem- 
perature T/Tc = 0.05 (full line), 0.25 (long-dashed line), and 0.50 (short-dashed 
line) . 

The density profiles obtained from Eqs. (l)-(4) for the gas at three differ- 
ent values of the temperature (in units of the critical temperature = 
{\/6N/nkB)hw± of the ideal Bose gas) are shown in Figs. 1 and 2. In the 
absence of a vortex (Fig. 1), the main point to notice is that the growth of the 
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Fig. 2. Density profiles for the condensate and for the thermal cloud in the presence 
of a vortex. Units and symbols are as in Fig. 1. The inset shows an enlarged view 
of the profiles near the center of the trap. 

thermal cloud exerts an increasing repulsion on the outer parts of the conden- 
sate, constricting it towards the central region of the trap. This effect is seen 
only when the condensate is treated by the NLSE and persists in the presence 
of a vortex (Fig. 2) but , as we shall see below, is missed in the Thomas- Fermi 
approximation where the radial kinetic energy of the condensate is neglected. 
In addition, the thermal cloud penetrates the core of the vortex and enhances 
the expulsion of the condensate from the core region (inset of Fig. 2), in a 
manner which again is governed in its details by the radial kinetic energy 
term in the NLSE. 
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Fig. 3. Density profiles for a BEC containing a vortex at T/Tc = 0.50 in the full 
calculation using the NLSE (full lines) and in the Thomas-Fermi approximation 
(dashed lines). The units are as in Fig. 1. The inset shows an enlarged view of the 
profiles near the center of the trap. 

The profiles obtained for a BEC containing a vortex at T = 0.5 by the 
full numerical calculation using the NLSE are compared in Fig. 3 with those 
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obtained in the Thomas-Fcrmi approximation. The constriction of the conden- 
sate in its outer parts and its expulsion from the core region by the thermal 
cloud are clearly underestimated in the Thomas-Fermi theory. 
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Fig. 4. The energy of a vortex, expressed as the thermodynamic critical frequency 
Qc in units of the radial trap frequency a;^, as a function of T/Tc. The results from 
the full calculation using the NLSE (full line) are compared with those obtained 
from Eq. (6) with the corresponding values of /i(T) (dotted line) and with the 
Thomas-Fermi values of /x(r) (short-dashed line). The long-dashed line shows the 
results obtained in a calculation using the Thomas-Fermi theory. 

Finally, Fig. 4 reports our results for the energy of the vortex as a function 
of T/Tc. The inaccuracies arising in the density profiles from the Thomas- 
Fermi treatment of the interplay between the condensate and the thermal 
cloud clearly lead to large errors in the estimation of flc{T). 

On the other hand the simple expression given in Eq. (6), with the two al- 
ternatives of using in it the chemical potential //(T") from the Thomas-Fermi 
approximation or from the full calculation using the NLSE, gives a reasonable 
account of the vortex energy up to T ~ 0.5 Tc- 



4 Summciry and future directions 

In summary, we have calculated the density profiles and the thermodynamic 
critical frequency for vortex nucleation in a strictly 2D Bose- Einstein-condensed 
gas at various temperatures. Our calculations have demonstrated the interplay 
between the thermal cloud and the structure of the condensate in the regions 
where the radial kinetic energy of the condensate is playing an important role 
(namely, in the outer parts of the condensate and in the vortex core) and 
validated the use of a simple expression of the vortex energy for temperatures 
up to about 0.5 Tc- 
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As to future directions, a practical realization of the 2D regime that we have 
investigated would require a combination of axial squeezing of the trap and 
an enhancement of the scattering length. In relation to the latter, a treatment 
of a Feshbach resonance in a 2D gas will be reported shortly. 
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